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ABSTRACT 

Studies of the broad emission-line regions (BLRs) in quasars have revealed solar or 
higher enrichment levels up to the highest redshifts. In combination with the presence 
of large amounts of dust in QSOs at z ^ 6, this implies that substantial amounts of 
star formation and nucleosynthesis took place at significantly earlier epochs. Here, we 
examine whether a top-heavy stellar initial mass function (IMF) is indicated by current 
data, by modelling the contributions from different regions of the IMF, including Type 
la/II and pair instability supernovae, to the metal synthesis in BLRs. We find that, 
in order to reproduce the observations of roughly solar values of N/C and Fe/Mg 
in these objects, (i) stars with a present-day IMF are sufficient, regardless of their 
metallicity, (ii) zero-metallicity stars with a top-heavy IMF severely underproduce 
N/C, and (iii) the contribution of Type la SNe is not strongly required by the data. 
Therefore, stars of mass ~ 1-40 Mq must have existed at z ~ 10-20, possibly coeval 
with any hypothesized stars of masses <; 100 Mq at these epochs. This is in agreement 
with the nucleosynthetic abundance pattern detected in extremely metal-poor stars 
in the galactic halo. 

Key words: cosmology: theory - first stars - galaxies: abundances - galaxies: starburst 
- quasars: emission lines 



1 INTRODUCTION 

Multiwavelength observations of increasing quality over the 
last decade have begun to strongly constrain the link be- 
tween the earliest epochs of star formation in the universe 
and the metal enrichment of high-redshift {z) protogalaxies 
and the intergalactic medium (IGM). Constraints on abun- 
dances at the highest redshifts currently comes from data 
on QSOs, where the chemical composition of the gas can be 
estimated from spectral features from the broad emission- 
line regions (BLRs) as well as broad and narrow absorption 
lines, probing regions ranging from sub-pc scales to those 
exceeding several tens of parsecs. 

Considerable progress has been made in recent years 
in observing and modelling such QSO environments, and 
detailed chemical evolution models have been developed to 
better assess reliable abundance diagnostics (Hamann & Fer- 
land 1999). The NV/CIV line ratio has received particular 
attention because of its high sensitivity to the overall metal- 
Ucity of the gas. For typical values of BLR parameters, the 
models predict NV/CIV ~ 0.1 for solar N/C. The current 



data indicates that most QSOs have supersolar N/C and 
BLR gas metallicities, including the highest-^ ones at 2; ~ 
6.28 and 5.99, which have NV/CIV = 0.35 and 0.67 respec- 
tively (Pentericci et al. 2002). In addition, a roughly con- 
stant value of Fell/Mgll (~ 2-15) is detected in QSOs up to 
the highest redshifts at z ~ 6.4 (Freudling, Corbin & Korista 
2003; MaioUno et al. 2003; Dietrich et al. 2003a and refer- 
ences therein). Assuming that Fell/Mgll traces Fe/Mg and 
using Fell/Mgll ~ 3 for solar Fe/Mg, the data imply BLR 
gas metallicities of near- or exceeding solar values (Dietrich 
et al. 2003b). Finally, recent submm/mm/IR observations 
of high-z QSOs have revealed the presence of large amounts 
of dust (thermal dust masses ~ IO^Mq) up to redshifts as 
high as 6.4 (Bertoldi et al. 2003), implying associated star 
formation rates of ~ 10^ M© yr"'^. 

Although the sizes and masses of the physical regions 
probed by BLR and dust observations are orders of magni- 
tude apart, together they yield a consistent scenario: that 
significant stellar activity and nucleosynthesis occurred at 
2: ^ 6. At these epochs, the age of the universe is approx- 
imately 1 Gyr and approaches the enrichment timescale of 
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Table 1. Properties of the three stcUar IMF models considered 
here: Miow and M^igh arc the lower and upper mass limits in 
Mq of the stellar mass range, a is the exponent of the power- 
law IMF (see text), is the initial stellar metallicity and the 
last entry indicates whether the contribution from Type la SNe 
is considered. 



Model 


Mow 


Mhigh 


a 


2* 


Type la 


A 


1 


100 


2.35 





no 


B 


1 


100 


2.35 




yes 


C 


100 


1000 


1 





no 



low- mass stars and Type la supernovae (SNe). These data 
may imply that star formation in such environments oc- 
curred rapidly and involved an initial mass function (IMF) 
favoring massive stars (Hamann & Fcrland 1993). Our aim 
here is to test whether the observed abundance ratios in 
combination with the constraints from stellar evolutionary 
timescales can place interesting limits on the nature of the 
stellar IMF at early epochs. In particular, the presence of 
significant quantities of iron has led some authors to pro- 
pose a Type la SN origin. Alternatively, Type II SNe and 
pair instability SNe (PISNe; Woosley & Weaver 1995; Heger 
& Woosley 2002) from stars of respective mass ranges 10-40 
M0 and 14O-16OM0 could also be substantial sources of iron 
in QSO BLRs. These massive stars have the added advan- 
tages of providing the same nucleosynthesis site for alpha- 
elements such as magnesium, and of being sources that can 
create dust and generate prompt metal enrichment. 

In this work, we examine the evolution of N/C and 
Fe/Mg for cases involving a top-heavy and present-day IMF 
in a burst scenario at early epochs. Our approach here is 
simple and chosen to reveal whether a preferred stellar mass 
range is indicated by the available data. For more detailed 
modelling of various star formation scenarios, we refer the 
reader to, e.g., Romano et al. (2002), Hamann & Ferland 
(1999) and references therein. 



2 STELLAR IMF MODELS 

We take the stellar IMF to have the form, 0(M) = (j)oM~", 
where a and M are the IMF slope and stellar mass respec- 
tively. We consider two cases: (1) a present-day form, where 
a equals the Salpeter value of 2.35, and the stellar mass 
range is 1-100 M©, and (2) a flat top-heavy IMF involv- 
ing very massive stars (VMSs) with a = 1 (ensuring a flat 
mass- weighted IMF, M(j}{M)) in the mass range 100-1000 
Mq . We have chosen this division in order to distinguish the 
stellar signatures of these two IMFs, although it is possible 
that these IMFs were coeval in the past, given the right 
combination of conditions (Scannapieco, Schneider & Fer- 
rara 2003). Both cases are normalized over their respective 
mass ranges as, / dMM0(M) = 1. For the first case, we 
consider two stellar metallicities, = and = Zq , and 
only = for the second, since it is currently thought 
that stars preferentially form in a top-heavy IMF with char- 
acteristic masses of a few 100 Mq when the gas metallicity 
is below a critical value of Za = 10~^^^Zq (Bromm et al. 
2001; Schneider et al. 2002, 2003a). 



The adopted ejected masses in C, N, Mg and Fe as a 
function of stellar mass and metallicity are collected in Ta^- 
ble 2. We follow Venkatesan & Truran (2003) and use van 
den Hoek & Groenewegen (1997) for stellar masses < 8 M0. 
Since no detailed compilation exists for Z^ = intermedi- 
ate mass star (IMS) yields, we approximate this by using 
= 0.001 (0.05 Zq) IMS yields for = 0, which do not 
appear to be substantially different by mass fraction (see 
Venkatesan & Truran 2003 and Chieffi et al. 2002). The 
yields from Type II SNe for = Zq (Z* = 0) stars in the 
ll(12)-40 Mq mass range are taken from the case A mod- 
els from Woosley & Weaver (1995), corresponding to typical 
SN explosion energies of ~ 1.2 x 10'''"'^ erg. We assume that 
stars of mass 40 100 Mr.) implode directly into black holes 
and therefore neglect their metal yields. For our top-heavy 
IMF case, only the stars of mass 140-260 Mq avoid com- 
plete collapse into a black hole (Heger & Woosley 2002) and 
contribute to the nucleosynthetic output. We note that our 
adopted yields do not include the efl'ects of stellar rotation 
which could potentially be important, particularly at low 
metallicities (Meynet & Maeder 2002; Chiappini, Matteucci 
and Meynet 2003). Our intent here is to consistently com- 
pare the yields as a function of stellar mass and metallicity. 

For Type la SNe, we use the revised W7 models from 
Nomoto et al. (1997), which are essentially independent of 
the progenitor mass. We assume that 5% of the stellar IMF 
in the mass range 3-16 Mq is in binary systems which lead to 
Type la SNe through the single degenerate scenario (Gibson 
1997; Matteucci & Recchi 2001), and that the mass ratio 
distribution of the binary's stars is close to unity, with the 
metal input from Type la SNe occurring on the evolutionary 
timescales of the secondary star. Recent theoretical studies 
(e.g., Kobayashi et al. 1998) indicate that low-Z ([Fe/H] 
-1) progenitors experience an inhibition of the Type la SN 
phenomenon. Thus, we assume that Type la SNe do not 
occur for = stars. In combination with the two IMF 
cases above, this leads to the three models (A, B and C) that 
we consider in this paper, as detailed in Table 1. Although 
there remains significant uncertainty in the exact values of 
stellar yields as a function of mass and metallicity, these 
cases should roughly bound the contributions from various 
stellar sites of origin. 

For the elements considered here, we include the total 
abundances of all the relevant isotopes of each element in 
its final yield; for iron, this additionally includes ^®Ni and 
^^Ni which eventually decay into iron isotopes. Although the 

isotopes of most elements are rare compared to the dominant 
form, some exceptions of interest here are N, Mg and Fe. We 
additionally note that the term yield in this paper connotes 
the total ejected mass fraction in individual elements or in 
metals; in the case of IMSs, since the yields from van den 
Hoek & Groenewegen (1997) correspond to the net metal 
output, we account for the original metal composition of 
the star as well. 

Other than the input from Type la SNe which is mod- 
elled as described above, the timescales for stellar yields cor- 
respond directly to the main-sequence stellar lifetimes for a 
burst mode of star formation. That is, for a single burst at 
z = 30, the yield at 2; = 6 equals the cumulative yield from 
those stars whose lifetimes do not exceed the time difference 
between z = 30 and z = 6. The main-sequence lifetimes are 
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Figure 1. Ejected mass fractions Yi = Mi/M as a function of 
stellar mass, M, for two values of Z*, where i can be C, N, Mg 
and Fe. Dashed lines (circle points) and solid lines (square points) 
represent Yi for Z* = and Z* = Zq stars, respectively. 



taken from Schaller et al. (1992) for Zi, = Zq stars and from 
Schaerer (2002) for Z = stars. 



3 RESULTS 

Figure 1 displays the ejected mass fractions Yi = Mi /M as a 
function of M for two values of (0 and Zq), where i can 
be C, N, Mg and Fe. We see that carbon has contributions 
from the entire IMF, whereas the production of N is primar- 
ily limited to the part of the IMS range (^ 3-4 Mq) that 
corresponds to efficient horizontal branch burning (Henry 
2003). Mg and Fe production, not counting Type la SNe, is 
limited to stellar masses ^ 10 Mq : Type II SNe eject compa- 
rable mass fractions in Mg and Fe, whereas PISNe produce 
an almost mass-independent iMg and a Vpe that rapidly in- 
creases with stellar mass. Figure 1 also shows the well-known 
trend that metal production increases with decreasing for 
IMSs (Henry 2003) but rises with increasing for massive 
stars (Woosley & Weaver 1995). For comparison, the yields 
available for Z-, = IMSs (Chieffi et al. 2002) for Reimers 
mass loss parameter r) = Z are Yc = [5.3,3] x 10"'' and Fn 
= [3,4] X 10~^ for stellar masses 4 and 7 Mq respectively. 

In order to facilitate direct comparison with observa- 
tions. Figures 2 and 3 show the evolution of the element 
ratios Fe/Mg and N/C by mass as a function of redshift and 
cosmological age for the three models A, B and C specified 
in section 2. Also shown in each figure is a shaded rectan- 
gle indicating the approximate observed range of BLR gas 
metallicities from studies of each abundance ratio, ^ 1-10 
Zq from NV/CIV over z ~ 0.0-6.3 (Hamann & Ferland 
1999; Pentericci et al. 2002), and ~ 1-5 Zq from Fell/Mgll 
over z ~ 0.1-6.4 (Freudling et al. 2003; Maiolino et al. 2003). 
The solar values for N, Fe, and Mg were taken from the re- 
vised photospheric abundances in Holweger (2001) and that 
for C from AUende Prieto, Lambert & Asplund (2002), so 



that (Fe/Mg)0 ~ 1.9 and (N/C)© 0.4 by mass. The rect- 
angle in the plots is intended to provide a rough reference for 
comparison between the models here and the current data 
indications of ~ solar or higher values of N/C and Fe/Mg 
in QSO BLRs. We emphasize that this box is not an aver- 
age of the data, which show significant scatter at a given 
epoch and which measure the ratios of the observed fiuxes 
in specific element ionization states. Converting these, e.g., 
Fell/Mgll, to the underlying abundance ratios (like Fe/Mg) 
involves the modelling of complex BLR environment physics 
(Hamann & Ferland 1999) and related uncertainties, which 
we have not accounted for here. 

For models A and B, each set of curves for a specific Zi, 
corresponds to the burst mode of star formation for three 
burst turn-on redshifts of 30, 20 and 10. These are chosen 
to be consistent with the range of formation epochs of first 
stars from the W MAP data on the microwave background 
(Spergel et al. 2003) and from cosmological simulations (Cia- 
rdi, Ferrara & White 2003). We follow the evolution for as 
long as it takes all the stars in each IMF to evolve after the 
burst. All the curves in Figure 2, and that corresponding to 
model C in Figure 3, are artificially extended to z = for 
comparison at z 6. The yields from VMSs are assumed to 
be instantaneously generated, given the average VMS life- 
time of ~ 2 Myr {Az ~ 0.4 at z ~ 30). 

Figures 2 and 3 show that massive stars and VMSs con- 
tribute trace amounts of elements like N at high redshifts. 
As noted above, N, in contrast to C, is primarily made from 
a narrower mass range in metal-poor or metal-free stars, 
about 4-7 Mq. In Figure 3, we see that for models A and B, 
after some initial fiuctuations owing to the detailed yields 
at the upper end of the IMF, N production peaks when the 
burst age is about 10® yr, corresponding to the lifetime of 
a 4 Mq star, the longest timescale on which N is created 
in this case. C production, on the other hand, continues to 
rise steadily, leading to the eventual steady decline of N/C. 
For Fe/Mg in Figure 2, the situation is simpler for models A 
and C where the prompt element generation from ^ 10 Mq 
stars is seen. The curves for model B in Figure 2 display 
the increasing contribution of Type la SNe at late times, 
peaking at cosmological ages corresponding to the time lag 
for the evolutionary timescale oi a, Z — Zq 1.5 Mq star, 
the smallest mass for the secondary star in the Type la SN 
model that we consider here. 

Our results may be summarized as follows: (i) model 
C (a top-heavy IMF with PISNe) can explain the observed 
Fe/Mg values but not N/C, due to the severe underproduc- 
tion of N by PISNe. Note that the latter conclusion may 
change with future revisions of VMS models; as Heger & 
Woosley (2002) point out, their exclusive study of helium 
cores may not adequately track N production by VMSs. (ii) 
Models A and B, which represent a present-day IMF with 
nucleosynthetic contributions from Type II and/or Type la 
SNe, can reproduce the observed ranges of both Fe/Mg and 
N/C at near- or supersolar values, regardless of stellar metal- 
licity. (iii) The inclusion of Type la SNe pushes the predicted 
Fe/Mg towards the upper end of the observed values (each 
Type la SN event generates Fe/Mg ~ 87 and N/C ~ 2.4 
xlO~^), but is still within the observed range of Fe/Mg, 
as shown by some earlier papers (see section 1). (iv) Even 
if alternate Type la SN mechanisms are considered that 
lengthen the creation timescale beyond 1 Gyr (see section 
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Figure 2. Evolution of Fe/Mg as a function of redsliift and cos- 
mological age. Dashed, solid and daslied-dotted lines represent 
models A, B and C respectively (see Table 1). The three curves 
in models A and B represent burst turn-on redshifts of z = 30, 
20 and 10. The shaded box indicates the approximate observed 
range of BLR gas metallicities of 1-5 Zq from Fell/Mgll over 
2 ~ 0.1-6.4. 
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Figure 3. Same as Figure 2, but for N/C. The shaded box indi- 
cates the approximate observed range of BLR gas metallicities of 
1-10 Zq from NV/CIV over z ~ 0.0-6.3. 

2), models A and B would still successfully reproduce the 
data by 2 ~ 6, i.e., Type la SNe are not strongly required 
by the data. 

These conclusions should be fairly robust since we have 
considered the evolution of abundance ratios, which are in- 
dependent of parameters like the BLR mass, star formation 
efflciency, and the relative escape fractions of metals from 
the BLR versus fallback into the central black hole. 



4 DISCUSSION 

The study of element evolution in the deep potential wells 
associated with the BLRs near the most massive detected 



high-2 objects, rather than in the IGM (e.g., Ellison et al. 
2000), has the potential advantage of not being subject to 
the uncertainties of metal transport mechanisms. BLR met- 
als may therefore reflect the original element synthesis more 
accurately, barring any enriched gas that falls into the cen- 
tral black hole. Both IGM and BLR studies share, however, 
sizeable modelling uncertainties in translating observed line 
strengths to the underlying metallicity, which depends on 
the gas temperature, density and ionization balance. 

Within these uncertainties, our results indicate that 
stars in a present-day IMF, regardless of stellar metallicity, 
can reproduce the data on those QSO BLR abundance ra- 
tios that we have considered here. A new population of stars 
(top-heavy IMF or otherwise) is not required from such data, 
although VMSs at z ~ 10-20 may be found to be attractive 
candidates to explain other observations, e.g., in order to 
match the observed amplitude and anisotropy of the near- 
infrared background (Salvaterra & Ferrara 2003; Maglioc- 
chetti, Salvaterra & Ferrara 2003; Cooray et al. 2003) and/or 
if the high Thomson optical depth initially detected by 
W MAP is not revised by future data. 

Therefore, the epochs at z ^ 6 may not have been ex- 
clusively dominated by a top-heavy IMF and in fact appears 
to require the presence of stars in a present-day IMF in addi- 
tion to any other coeval stellar population. Indeed, a prompt 
transition from a top-heavy IMF to a Salpeter IMF may well 
have occurred in the dense environments probed by QSO 
BLRs, as the metals released by the first episodes of star 
formation could remain confined and efficiently pollute the 
gas up to Zcr. Not requiring an exclusive VMS contribution 
by 2 '-^ 6 in QSO BLRs would also be consistent with, e.g., 
the mild evolution of Fe/Mg over z ~ 0.1-6.4, since metal- 
free star formation is not observed at these lower redshifts. 

If Type la SNe are assumed to be the dominant sources 
of iron, their generation timescales in relation to the BLR Fe 
abundances at 2; ~ 6 places the first star formation epochs at 
z - 10-20 (Freudling et al. 2003), consistent with WMAP's 
findings. However, it is unclear if Type la SNe can occur 
when Zi, = Q stars are involved, and consequently, barring 
extremely short metal enrichment timescales, whether such 
SNe play a significant role in early cosmic metal generation. 
Observations of large amounts of dust in the highest-z QSOs 
(Bertoldi et al. 2003) provide another avenue to constrain 
early SNe. Although the connection between the synthesis 
timescales of metals and dust is not yet completely under- 
stood, it is currently believed that dust could be created 
from Type II SNe (Todini & Ferrara 2001) and/or PISNe 
(Schneider, Ferrara & Salvaterra 2003; Nozawa et al. 2003), 
but not from Type la SNe. As Bertoldi et al. (2003) addi- 
tionally point out, there may not be enough time to create 
dust from intermediate-mass stars by z ~ 6. 

Finally, the ashes of the earliest generations of stars 
may be reflected in the composition of very metal-poor halo 
stars in the local universe. Observations of such low-mass 
stellar relics in our Galaxy reveal the presence of elements 
heavier than Mg which cannot be synthesized in their inte- 
riors, and which must therefore have been already present 
in the pre-enriched gas cloud from which such stars formed. 
The recent detection of the most iron-poor star observed to 
date (Christlieb et al. 2003) lies tantalizingly close to the 
hypothesized Zcr for the formation of VMSs (see section 
2), although this star is relatively enriched in other metals 
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such as C. This star's elemental abundance pattern indicates 
that the original starforming cloud was enriched with metals 
ejected by SNe from ~ 20-130 Mq stars (Umeda & Nomoto 
2003, Tumlinson, Venkatesan & ShuU 2004), and it is only 
partially consistent with enrichment by PISNe (Schneider 
et al. 2003a). This derived mass range for a first-generation 
stellar population is roughly similar to that which could be 
responsible for the enrichment of QSO BLRs as shown in 
this work. Thus, a top-heavy IMF may not be required to 
explain the metal abundance patterns in QSO BLRs, which 
is consistent with the current data on extremely metal-poor 
stars in the galactic halo. It remains to be seen whether 
future observations of such local stellar relics and of the 
highest-z objects yield a consistent picture of the sources of 
the earliest metal synthesis in the universe. 
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